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Abstract 

This quarterly report deals with the development of a two-dimensional 
analytical model that describes the dynamics of an n-mass vertical tethered 
system. Two different approaches are described: in the first one the control 

quantities are the independent variables while in the second one the Cartesian 
coordinates of each mass expressed in the orbiting reference frame are the 
independent variables. The latter model has been used in the 3-mass version to 
simulate the dynamics of the tethered system in applications involving the 
displacement of the middle mass along the tether. In particular the issues 
related to reproducing predetermined acceleration profiles and "g-tuning" are 
reported . 


PRECEDING PAGE BLANK NOT FILMED 
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Figure Captions 


Figure 2.1.1 Geometry and generalized coordinates for the n-mass 

tethered system . 

Figure 2.3.1 Geometry and Cartesian coordinates for the n-mass 

tethered system . 

Figure 2.3.2 Schematic of the tension forces acting on a single 

mass of an n-mass tethered system. 

Figure 2.3.3 Geometrical relations between generalized coordinates 

and Cartesian coordinates . 

Figure 2.5.1 a-1 Dynamic response of the 3-mass tethered system for a 

sinusoidal motion of the middle mass with a peak-to-peak 
amplitude of 2000 m and an angular frequency of 10' 3 
rad/sec . 

Figure 2.5.2 a-1 Same as in Figure 2.5.1 with a peak-to-peak amplitude of 

the motion of the middle mass equal to 8000 m. 

Figure 2.6.1 a-m Dynamic response of the 3-mass tethered system during a 

transfer maneuver of the middle mass designed for achieving 
a predetermined, steady state g-level ("g-tuning") . The 
control law is a half cycle sinusoid with a peak-to-peak 
amplitude of 2000 m. 

Figure 2.6.2 a-m Same as in Figure 2.6.1 with a control law of different type. 

A modified hyperbolic tangent is used to perform the transfer 
maneuver instead of the half cycle sinusoid. 
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1.0 INTRODUCTION 

This is the third Quarterly Report submitted by SAO under contract NAS8- 
36606, "Analytical Investigation of the Dynamics of Tethered Constellations in 
Earth Orbit (Phase II)," Dr. Enrico Lorenzini, PI. It covers the period from 22 
September 1985 through 21 December 1985. 


2.0 TECHNICAL ACTIVITY DURING REPORTING PERIOD AND PROGRAM STATUS 

2.1 Two-Dimensional Equations Of Motion Eor The N-Mass System 

The Lagrangian formulation has been used to derive the two-dimensional 
equations of motion for an n-mass tethered system. The equations of motion are 
expressed in the orbiting reference frame formed by the local vertical (z) and 
the local horizon (x) centered at the system C.M. The Lagrangian coordinates, 
shown in Figure 2.1.1, are the tether lengths (fj.) , the deflections (ej.) of each 
mass from the line connecting the first mass to the last mass and the angle ( 6 ) 
between that line and the local vertical . The assumptions are like those used 
for the 3-mass system, namely: point masses, circular orbit of the C.M. , second 

order gravity potential expansion, elastic tethers. Following a procedure simi- 
lar to the one in Quarterly Report #1 and referring to Figure 2.1.1, the kinetic 
energy of the system is given by: 


i n 

T - IE 


mil Vji | 2 


( 2 . 1 . 1 ) 


In equation (2.1.1) the velocity of the i th -mass is given by: 
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vii = v ± + (Q x r i) 


and since we have : 


v± = v Bi + (6 x p Bi ) 


r ± = a + pi 


we can write equation (2.1.2) as follows: 

Vii = v B1 + (fl x a) + (0+0) x p B1 


( 2 . 1 . 2 ) 


(2.1.3) 


(2.1.4) 


(2.1.5) 


so that: 


| Vii | = [x B i + Oa cos 8 + (0 — 0)z Bi ] -|- [z Bi + fia sin 9 — [8 — fl) x B i J (2 . 1 . 6) 


Since the body reference frame is centered at the system C.M. we have: 


n n 

mi x Bi = nii z B1 = 0 (2.1.7) 

i l 


Moreover, since the C.M. remains, through time, at the center of the orbiting 
reference frame we also have: 


n n 

Y m i x B i = Y. m t z B j = 0 
1 1 


( 2 . 1 . 8 ) 


By manipulating equation (2.1.6) we finally get: 
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T 


f mi {* B i 2 + ^Bi 2 
Z 1 1 

+ 2($ — fi) (xbiZbi — 


+ (0 -0) a (XB1 2 + ZBi*) + 

XbiZbi)} + J m tot n a a 2 


where the last term is the kinetic energy of the system concentrated 


The potential energy of the system is given by: 


where : 


n 



r± | = [xi 2 + (a - Zi) 2 ] V2 


We can therefore write the potential energy as: 


V 


m A [i + ^ - 


Xbi 2 + z Bi 2 3 Zj 2 1 

2a 2 2 a 2 J 


where : 


xbi 2 + Zbi 2 


2 2 
Xi + Zi 


Since : 


H = tl 2 a 2 


zi = z B i cos 9 - x Bi sin 6 


(2.1.9) 
at the C.M. 

(2.1.10) 

( 2 . 1 . 11 ) 

( 2 . 1 . 12 ) 

(2.1.13) 

(2.1.14) 


after algebraic manipulations and because of equations (2.1.7) we finally get: 
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V = - | n>i{n 2 (3 cos 2 0 - l)z Bi 2 + n 2 (3 sin 2 0 - l)x B1 2 

-6 fi 2 sin 6 cos 6 x Bi z B1 )j - (2.1.15) 

In equation (2.1.15) the last term is the potential energy of the system when 
the total mass is concentrated at the C.M.; its value is constant for a circu- 
lar orbit. 

Assuming that €i « ti and defining Ri = mi/m to t we can write: 

z B i = “ Zbcmi 

z B j — ^i"Zbcmi 


i 

z B i ~ zbcmi 

J=i 


and : 


n-1 

ZBn = yi f-i ~ Zbcmi 
J=i 


Xbi — “ *BCM1 

XbJ = - X B CM1 


(2.1.16) 


(2.1.17) 

Xbi = € i Xbcmi 


x Bn — x B cmi 
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where zbcmi and x bcmi are the coordinates, in a body reference frame with the 
origin at mi, of the system C.M. They can be expressed as follows: 

zbcmi = R(a,...,n) ^1 + .... + R n ^n-i 


(2.1.18) 

n-2 

x bcmi Rj+i 

J=i 


where : 


R(i n) = (2.1.19) 

1=1 

Equations (2.1.16) and (2.1.17) relate the coordinates of the masses in body 
axis to the Lagrangian coordinates. 

The Lagrangian function (L) is readily obtained by subtracting the potential 
energy from the kinetic energy so that the equations of motion are given by: 



After many algebraic manipulations the final structure of the equations of mo- 


tion is as follows: 
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A[£ + bl] + 2(4 -Q). B e + (4 - d) C e = Q /m tot 


D[« + ge] + 2(4-0) E l + (0 + d) F t = Q /m tot (2.1.21) 


G[4, 2(4—0), 30 2 sin 0 cos 8 , 3fi 2 (sin 2 0 — cos 2 0 ) ] + eH t + the = Qe/m to t 

In equations (2.1.21) t, t , £ are column vectors of dimension n-1 and e, «, e 

are also column vectors of dimension n-1 with the last row equal to zero. The 
matrices A, B, C, D, E, E, H, L are non dimensional matrices and G is on the 

contrary dimensional, being a function of e, i, t, i. Their complex expressions 

are not reported in this quarterly report for the sake of brevity. The coeffi- 
cients in equations (2.1.21) are given by: 

b = fl 2 (1-3 cos 2 $ ) - (4-fl) 2 

d = 3n 2 sin0cos0 (2.1.22) 

g = n 2 (1-3 sin 2 5) - (4-0) 2 

The equations of motion (2.1.21) together with (2.1.22) will be fully developed 
in the next section for the case n=3 . Note that equations (2.1.22) are a set of 
2 (n-1) equations while in general the degrees of freedom of n masses in a two- 
dimensional space are 2n . In our case, however, the motion of the system C.M. 
is constrained so that the total degrees of freedom of the system are 2n-2. 
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2.2 Revised Two-Dimensional Equations Of Motion For The 3-Mass System 

From the generalized case of the previous section we derive in this section 
the equations of motion for a 3-mass system. The tedious computation of the 
matrix components in equation (2.1.21) is not reported. The equations of motion 
finally result in: 

Ri(l-Ri) [£i + b£ x ] + R X R 3 [£ a + b £ a ] - 2 (5-n)R 1 R 2 £- (£-d)R!R 2 e = Q £i /m tot 
R 3 (l-R 3 )[£ 2 + b£ a ] + R X R 3 [£ x + b£ x ] + 2(^-n)R 3 R 3 e +(^-d)R 3 R 3 c = Q /m tot 

c a 

R a (l-R a ) [ € + gc] - 2(^-n)R 3 v - (0+d)R 2 v = Q £ /m to t (2.1.23) 

£{Ri*i[*i + v] + R 3 £ 2 [£ 2 -v] +R 2 (1-Rj) 6 2 j + 

+ 2(d-0){R x £i.[ii+v]+Rj£ 2 [£ 2 -v]+R 2 (l-R 2 )«e} + 

+ 30 2 sinflcose|R 1 £ 1 [£ x +v] +R 3 £ 2 [£ 2 -v] -R a (1-R a ) e 2 } + 

+ 3D 2 (sin 2 0-cos 2 0)R 2 ve-R a £{R x £ x -R 3 £ 2 }-R 2 ev = Qg/m t0 1 

where b, d and g are given by equations (2.1.22) while v and v are given by: 

V = R-3^2 “ Rl^l 

v = R 3 ^ a - R x i x (2.1.24) 

The generalized forces had already been derived in Quarterly Report #1; they 
are given by: 

°£x = 

Q = -T 3 (2.1.25) 

«• 2 

Q« = ~c(Ti/£x + T a /£ a ) 


Qe 


0 
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A comparison of equations (2.1.23) with equations (2.1.15) of Quarterly Report 
#1 , taking into account the different notation, shows that some inaccuracies 
were present in the latter equations. However those inaccuracies involved terms 
which multiplied v or v . v and v represent the distance of m a from the system 
C.M. and its rate of change respectively. In all the simulations reported in 
Quarterly Reports #1 and #2 the values of v and v were zero or close to it so 
that the results shown in those reports are correct. As a confirmation of this 
we ran with the new model one of the case reported in Quarterly Report #2 
obtaining the same results. When simulating the middle mass travelling along 
the tether the incorrect terms are, on the contrary, important and the use of 
the new model is therefore mandatory. 

In equations (2.1.23) the independent variables are the controlled or the 
observed parameters. To associate control laws to the system equations is 
therefore a straight forward process. Equations (2.1.23) , however, have the 
disadvantage of being non-normalized . In matrix form the equations (2.1.23) can 
be written: 


A x = b (2.1.26) 

where x = [ 6 , c, l lt £ a ] T , A is a symmetric 4x4 matrix and b = [b a , b a , b 3 , 
b,»] T . The numerical integration of equations (2.1.26) implies the inversion of 
matrix A. Even if it is not a problem for a 4x4 matrix, the inversion becomes 
more and more time consuming for bigger matrices (greater number of masses) . 
For this reason a different model in Cartesian coordinates has been derived, as 
shown in the next section. Eor the sake of completeness the components of 
matrices A and b are given below: 
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an 

— 

R x £i 

[£ x +v] +R3C2 [^a-v] +R a (1-R a ) e 

a 22 

= 

R a (l 

-R2) 

a 33 

= 

Ri(i-Ri) 

a 4 4 

= 

R 3 (1'Rj) 

a 12 

= 

a 21 

= -R 3 v 

a 13 

= 

a 31 

= -R x R a e 

a 14 

= 

a 4i 

= RjR 3 6 

a a a 

= 

a 32 

= 0 

a 3 4 

= 

a 42 

= 0 

a 34 

= 

a 43 

= R1R3 


-2 (4-fi) {r^ [i 1+ v] +R 3 £j [£ 3 -v] -t-Rj (1-Rj) ee}- 

d|Ri^i[^i+v] +R 3 £ 3 [£ 3 -v] -Rj( 1-R 3 ) e 2 }-30 2 (sin 2 0-cos 2 0)R x ev 

Qe/mtot " R a (l-R a )g e + 2(^-0)R 3 v + dR a v (2.1.28) 

Qa /nitot ” b [R x (1 - R X ) £ x + R X R 3 £ 3 ] + 2($*"0)R X R 3 £ - dRiR 3 € 

Q, /mtot - b[R 3 (l-R 3 ) £ 3 +R x R 3 £ x ] - 2(0-D)R 3 R 3 e + dR 3 R 3 e 
<-a 

Equations (2.1.27) and (2.1.28) are complemented by equations (2.1.25), (2.1.24) 


bi = 

b 3 = 
b 3 = 
b 4 = 


and (2.1.22). 
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2.3 Another Model For The N-Mass System. Cartesian Coordinate Equations Of 

Motion 

The model derived in the previous section is very convenient from the con- 
trol viewpoint since the variables of the system equations are the control pa- 
rameters. However, because of the difficulty of inverting large matrices, this 
model's numerical integration is not efficient. In this section, we will there- 
fore develop the equations of motion for a two-dimensional , n-mass tethered 
system in Cartesian coordinates relative to the orbiting reference frame. The 
equations in this form can be directly integrated without matrix inversion; on 
the other end the relations between the system variables and the con- 
trolled/observed variables are much more complex. 

Referring to Figure 2.3.1 let us define the radius vector of the generic 
i^-mass as follows: 


= a + pi. 


(2.3.1) 


where : 


Pi = xi I + Zi k 



(2.3.2) 


In the case of Q = constant we can write the inertial acceleration of the i**- 
mass as follows: 


?n = pi + 2fj x pi + 0 x (Ox ?i) 


(2.3.3) 
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Figure 2.3.1 
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This equation can also be expressed as: 


?ii = pi + 20 x pi + n (n ?i) - ?i 


_ 2 

n 


(2.3.4) 


Since in the two dimensional case fi is always perpendicular to r A the third term 
in equation (2.3.4) is equal to zero. By using equations (2.3.2) we finally 
have : 


?n = [*i - 2flzi - 0»Xi]l + [zi + 2flxi — 0 2 (Zi-a)]k (2.3.5) 

The equation of motion for the i th -mass is therefore: 


nurii = F gl + F T i + F pi 


(2.3.6) 


In equation (2.3.6) F g i denotes the gravity force acting upon the i^-mass while 
Fti is the net force due to tether tension acting upon the same mass and F pl is 
the perturbation force . 


Without limiting the generality of the model we assume F p i = 0. The gravita- 
tional force is obtained from the potential energy as follows: 


F g i = - V Vi 


(2.3.7) 


where Vi can be readily derived from equation (2.1.12) obtaining: 


Vi = | mi[20 2 a 2 + 2fl 2 aZi - fl 2 (xi 2 + z A 2 ) + 3n 2 Zi 2 ] 


(2.3.8) 


Substituting in equation (2.3.7) and performing the derivatives we get: 
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F,i = mitc-flzxi)! + (0 2 a + 20 2 Zl )k] 


(2.3.9) 


The two scalar equations of motion for the i^-mass can be finally written as 
follows : 


Xi — 20 zi — Fixi/mi 


(2.3.10) 


z t - 30 2 zi + 20 xi = Fixi/mi 


These equations are usually called the Hill's equations. 


Referring to Figure 2.3.2 


the force F Ti is given by: 


Fxi = - [Txi + T x ,i-i]I - [T xi - T. tl .x]k 


(2.3.11) 


where T x i , T xl are the components of the tension Tj in the tether connecting 
the i^-mass to the i+l^-mass; they are respectively given by: 


Txi = T t cos ^ = Ti 


Xj+i - x t 


[(* 1+1 + Xi) 2 + (Zi^i + Zi) 2 ] 


1/2 


(2.3.12) 


T x i = Ti sin£i = Tx 


2l + l ~ Zj 


[ (Xi^l - Xi) 2 + (Zi+i - Zi) 2 ] 


1/2 


Referring to Quarterly Report #1 the tension Ti is given by: 


Ti = ^ (£i - Z„i - to.) = k ti £ t 


(2.3.13) 


In equation (2.3.13) , £i is the geometric distance between the i^-mass and the 
i + l^-mass, £ ci is the controlled length of that tether segment, £di the length 
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of the associated longitudinal damper and £ t i the tether stretch. The length 
l C i can follow a prescribed control law if that particular tether segment is 
actually controlled or it is the tether length at rest if that tether segment is 
uncontrolled. The length £ di is obtained by solving the following equation: 


£di = ( £ i — £di — £ci) 

kdi 


kdi 

kdi 


£di 


(2.3.14) 


Equation (2.3.14) is the generalization of equation (2.2.5) in Quarterly Report 
#2. Kdi and K d i are the stiffness and damping coefficient of the i th - longitudi- 
nal damper respectively. 

We must now derive the expressions that relate the system quantities 

(xi x n , Zi z n , Xi,...x ni Zi i n ) to the output quantities ( 6 , 

£n-...£n-l. <i» l i...., £ n -i. £n-j) • In this report we will 

derive the expressions that relate the variables while the first derivatives, 
when necessary, are numerically computed in the computer program. Referring to 
Figure 2.3.3 we have: 



= [(Xi+l-Xi ) 2 + (Zi + 1 - Zi) 2 ] 1/2 


(2.3.15) 


The lateral displacement £i is derived by computing the coordinates of the point 
of intersection between the straight line through mi and m n and the straight 
line perpendicular to it passing through mi. The equations of the above men- 
tioned lines are respectively given by: 


z — Zi 


Z n ~ Zj 
Xn “ Xi 


(X - Xi) 


Xn ~ Xi 


z — z i+1 


Z n - Zi 


(X - Xi*i) 


(2.3.16) 
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The solution of equations (2.3.16) are the coordinates of the point of intersec- 
tion x c i, z ci . After some manipulations we get: 

X ci = [Xl - Xi^l t g 2 fl + (Z i + 1 - Zl)tg^)]/(1 - t g 2 ^) 

(2.3.17) 

Zci = [x c i - X 1 + 1 ] t g 6 + Z i+ 1 
So that finally is given by: 

«i = [(x 1+ i - x cl ) 2 + (z i+ i - z ci ) z ] V3 'sign (x i+i - x ci ) (2.3.18) 

Equation (2.3.18) and equations (2.3.17) provide the relations between system 
variables and output variables. 

2.4 The 3-Mass System. Cartesian Coordinate Equations Of Motion 

The equations derived in the previous section for an n-mass system are 
reported, in this section, for a 3-mass system. The equations of motion become 
as follows: 


x x = 2f) Zi + Ti (xj-xi) / (rai^i) 

x a = 2Q zj + Tj (x 3 -x 2 )/(m 2 £ 2 ) + Tx (xx-x 3 ) / (m 2 £x) 
x 3 = 2fiz 3 + T 2 (x 2 -x 3 )/(m 3 £ 2 ) 

z 2 = 3n 2 zx-2nxx + Tx (z 2 -zx) / (mx£x) (2.4.1) 

z 2 = 30 2 z 2 -2nx 2 + T 2 (z 3 -z 2 )/(m 2 £ 2 )+T 1 (zx-z 2 )/(m 2 £i) 
z 3 = 3fl 2 z 3 -2fix 3 + T 2 (z 2 -z 3 ) / (m 3 £ 2 ) 
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where : 


Ti = EiAi^ti/^ci 
Tj = EjAj£ t a/^ca 

£ t i = £i“^di"^ci (2.4.2) 

^ta = £a - ^da“^ca 

idl = J<tl£tl/l<dl _ l<dl^dl/l<dl 

ida = kta^ta/kda ” kda^da/kda 

The relations between system variables and output variables become: 

0 = tg'^ (xx-x 3 )/(zx-z 3 ) ] 

= [(xa-xx ) 2 + (za-zi) a ] 1/J 

l 2 = [(x 3 -x 2 ) 2 + (z 3 -z 2 ) 2 ] 1/J (2.4.3) 

«x = [ (xa~x q i) 2 + (z 3 -z cl ) 2 ] 1/a 'sign(x a -x cl ) 


where : 

Xci = [xi-x 3 tg 2 0 + (za-zx) tg0 ]/ (1- tg 2 0) 

z c i = [x c x-xa]tgfl + z 2 *-■- (2.4.4) 

Equations (2.4.1) through (2.4.4) have been numerically integrated in a newly 
developed computer code to obtain the results illustrated in the following sec- 


tions . 
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2.5 Generating Sinusoidal Acceleration Profiles 

A unique feature of a 3-mass, vertical, tethered system is the capability 
of providing a variable-g level in the middle platform by simply moving the 
platform along the tether. The middle platform experiences a vertical accelera- 
tion, linearly dependent upon the distance from the system C.M. , given in g-unit 
by the following formula: 


a v = 9.81 x 3n 2 £ c .M. (2.5.1) 

By controlling the distance £c.m. it is possible, in principle, to produce a 
desired vertical acceleration profile. Actually the scenario is more complex: 
any displacement of the middle mass from a rest condition generates a transient 
response of the tethered system. Namely, the motion of the middle mass induces 
longitudinal oscillations of the elastic tethers and its non-zero velocity in- 
teracts with the orbital rate producing a Coriolis force that deflects, away 
from the vertical, the acceleration measured at the middle platform. During the 
acceleration or deceleration phase, moreover, the non-zero acceleration of the 
platform is added to the gravity gradient acceleration. The acceleration pro- 
file obtained is therefore a result of many parameters, its departure from the 
nominal profile depending primarily upon the control laws adopted and the effec- 
tiveness of the dampers. 

The first case that we analyze is the generation of sinusoidal acceleration 
profiles with the middle mass starting from the system C.M. and moving in be- 
tween the upper end mass and the initial C.M. position. The Space Station mass 
adopted in the following simulations is consistent with the post-IOC phase of 
the Space Station program. The other parameters of the micro-g/variable-g labo- 



Page 24 


ratory are not intended to be definitive: they provide a possible configuration 
of the system that is not based upon optimization criteria. The system parame- 
ters are as follows: 


Space Station Mass (mi) : 

Middle Mass (m^) : 

End Mass (ms) : 

Tether Length : 

Orbital Height at the 
System C.M. : 


3.06752 x 10 5 kg 
5 x 10 3 kg 
10“ kg 
10 km 

450 km 


(2.5.2) 


The tether is a 2 . mm diameter kevlar tether. In these particular simulations 
the longitudinal dampers are non-adaptative ; they are tuned respectively to the 
longitudinal frequency of the associated tether segments at their initial tether 
length. The angular frequenices are given by: 


wtx — vWT^ml) 


Wl 2 = \/E A/ 


(2.5.3) 


where £ cl0 and £ c jo are the initial, controlled tether lengths for tether #1 and 
tether #2 respectively. The damping coefficient in both the longitudinal damp- 
ers is assumed equal to .9. The rotational (5) and lateral (e) active damping 
algorithms are operated according to the description given in section 2.3 and 
2.4 of Quarterly Report #2. 

The control laws adopted for the sinusoidal acceleration variation are as 


follows : 
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id = ^cio + A£ c I 1 + sin (w c t + 


i *>] 


-cJ 


= £ c jo - A£ c 1 + sin (w c t + 


I ”>] 


(2.5.4) 


In equations (2.5.4) A£ c is the amplitude of the sinusoidal displacement so that 
the peak-to-peak displacement is 2A£ C and w c is the angular frequency. Neglect- 
ing effects such as tether elasticity and lateral displacement of the mass m 2 , 
the resulting vertical acceleration measured at the middle platform, in zero 
order approximation, is given by: 


a v = 30*£ c>m . - A£ c Wc sin (w c t + - ir) 


(2.5.5) 


Assuming that the position of the C.M. does not move appreciably during the 
motion of the middle platform (as is the case in zero order approximation) we 
have: 


£c.m. s: A£ c 1 + sin (w c t + 


- *) 
2 


(2.5.6) 


Finally we get: 


a„ = 3n z A£ c + A£ c (3C1 2 — wf) sin (w c t + 


3 \ 
— ir) 

2 


(2.5.7) 


The resulting vertical acceleration variation is therefore dependent upon both Q 
and w c . The point of zero acceleration is no longer the system C.M. but occurs 
according to the following formula: 


w c t = sin' 1 [l/(w c /3n 2 - 1)] - | ir + k2ir 


k 


1 , 2 ,... 


(2.5.8) 
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If u e < O then the second term in equation (2.5.5) is negligible. Since the 
effects of tether elasticity and Coriolis acceleration are also smaller for low 
values of w c , the resulting acceleration variations of m a are closer to nominal. 

In the following set of simulations we assume w c = 10" 3 rad/sec while the 
effect of a different amplitude A i c is explored. The first group of plots shown 
in Eigure 2.5.1 a-1 is relevant to a peak-to-peak displacement 2A£ C = 2000 m. 
All the dampers are active. The simulation has been run for 11000 sec which is 
approximately 2 orbits. Eigure 2.5.1a shows the sinusoidal length variation of 
tether #1 vs. time. Eigure 2.5.1b represents the longitudinal damper length 
variation for the same tether. Figure 2.5.1c is the in-plane ( 6 ) angle vs. time 
while Figure 2. 5. Id is the phase plane 0-6. Figure 2.5.1e is the phase plane e- 
k where e is the lateral displacement of m a from the line through mi and m 3 . 
These last two figures show that the tethered system, after dissipating the 
initial transient motion, follows a steady state limit cycle sustained by the 
motion of m a . Figure 2. 5. If shows the components of the acceleration measured 
at the middle platform. The vertical component (a v ) is sinusoidal in shape but, 
as mentioned before, the zero-g acceleration point is no longer at the crossing 
of the system C.M. Because of the middle mass self acceleration the vertical 
component (a v ) overshoots towards the positive values (acceleration directed 
downwards). Note also that a much smaller, but not negligible, horizontal com- 
ponent (ah) appears. It is related to the Coriolis force. Figure 2.5.1g is an 
isometric, polar diagram of the acceleration. The vector connecting the 0,0 
point to any point along the curve gives the actual orientation and modulus of 
the acceleration. After the transient phase is over the acceleration vector 
rotates steadily around the C.M. of the middle platform at the angular frequency 
w c . Figure 2.5.1h is the isometric side view trajectory of the middle mass with 
respect to the system C.M. Negative numbers in the z axis meaning upwards. For 
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Figure 2.5. lgt 


Figure 2.5. lh 
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completeness Figures 2.5.1 i and 1 show the tension in tethers #1 and #2 respec- 
tively. 

A second case with a peak-to-peak tether length equal to 8000m is shown in 
Figure 2.5.2 a-1. Comments similar to those for Figure 2.5.1 hold true in this 
case. The shape of the trajectory side view of mass mj is less symmetrical with 
respect to the local vertical; in general, however, the behavior is qualita- 
tively very similar to the previous case. Quantitatively the acceleration com- 
ponents show a dependence with A£ e so that they are 4 times bigger than in the 
previous simulation . 


2 . 6 G-Tuning 

The possibility of varying the acceleration level in the middle platform 
provides a "g-tuning" capability. This means that the middle platform is moved 
from the zero-g point of the tethered system and is placed at the appropriate 
distance in order to have the requested g-level according to equation (2.5.1). 
Two different control laws are shown here. They are both open loop control laws 
in the sense that the acceleration in the middle mass is not fed back into the 
control system. The final acceleration value, therefore, depends only upon the 
final distance (£ c .m.) achieved. The first control law adopted is a sinusoidal 
control law, similar to the one illustrated in the previous section, that has 
been truncated after half a cycle. In Figure 2.6.1a-m the dynamic response of 
the system is shown in a case with a peak-to-peak displacement of the middle 
platform equal to 2000 m. Figure 2.6.1a depicts the tether length for tether #1 
vs. time: the sinusoidal control law is clearly visible. Figure 2.6.1b is the 
length variation of the longitudinal damper. Further comments on the figure 
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will be provided later on. Figures 2.6.1c and d are the in-plane angle vs. time 
and the phase plane d-Q respectively. These two figures show that after 7200 
sec there is still a small residual in-plane oscillation decreasing to zero. 
Similar comments are applicable to Figure 2.6.1e and f which represent the lat- 
eral deflection e and the phase plane e-e respectively. Figure 2.6.1g shows the 
acceleration components, along the orbiting reference frame, measured at the 
middle mass. In this "g-tuning" simulation the horizontal component disappears 
after the transient phase is over. The vertical component is however experienc- 
ing a steep variation at the time when the final distance from the C.M. has been 
reached because of the non-zero acceleration provided by the control law at that 
point. This effect accounts for the similar steep variation in the longitudinal 
damper's length. The system, therefore, needs additional time to reach a steady 
state acceleration condition. Figure 2.6.1h shows the isometric, polar diagram 
of the acceleration; it stresses again the negative effect of using a control 
law (sinusoidal) that has non-zero acceleration at the beginning and at the end. 
Figure 2.6.11 is an isometric plot of the trajectory side view of mass m 3 with 
the zero point located at the initial position of the system C.M. Figure 2.6.1& 
and m show the tension in tether #1 and tether #2 respectively. The larger 
oscillations in tether #2 tension are due to the expanded plotting scale. The 
results obtained by using the sinusoidal control law for "g-tuning" prompted us 
to adopt a control law with zero acceleration at the end of the maneuver. The 
new control law is a modified hyperbolic tangent; the tether lengths are con- 
trolled according to: 


£i = tio -I- A £ c (1 — e'^) tanh (at) 


( 2 . 6 . 1 ) 


£j = i -20 — A£ c (1 - e'fc) tanh (at) 



2400 


2200 


2000 


1800 


~ 1600 


O 

2 

X MOO 

f— 

u 

2 

U 

“? 1200 

a 

CxJ 

X 

E-* 

p 1000 


800 


600 


400 


X 

E- 

0 
2 
W 
-J 

01 

w 

CL 

2 

< 

e> 

d 

2 

O 


a 


1.95 


1.9 


1.05 


1.8 











THETAD (Deg/Sec) IN-PLANE ANGLE (Deg) 







EPSID (M/Sec) M3 LATERAL DEFLECTION, EPSI (M) 











M2 VERTICAL ACCEL..ACV (M/Sec**2) ( M2 ACCEL. COMPONENTS (M^/Sec*‘2) 









Tension Tether 1 (N) I l l l l , Z2 ,(M) 


200 


0 

-200 

-400 

-600 

-BOO 

1000 

1200 

1400 

1600 

1800 

2000 

2200 

400 

395 

390 

385 

380 

375 

370 

365 

360 

355 

350 










Page 45 


The two constants a and fi set the slope of the curve, /3 having been introduced 
in particular to slow down the acceleration phase (initial phase) of the control 
law. Results from a simulation with 2 A£ e = 2000 m, a = 1/2000 sec" 1 , = 1/1000 
sec" 1 are shown in Figure 2.6.2a-m. The sequence of plots is equivalent to the 
one previously illustrated. Note that the transient phase at the end of the 
maneuver has disappeared. In the same 7200 sec time frame the system reaches a 
quieter state than in the previous simulation. The initial transient phase is 
still there; it is partially due to the imperfect initial conditions. Figure 
2.6.2g shows the acceleration components in the local orbiting reference frame 
at the middle platform. Note that at the end of the maneuver the horizontal 
acceleration component has completely disappeared while the vertical component 
has reached its steady state value. 


2.7 Concluding Remarks 

The newly developed Cartesian coordinate model provides a very convenient 
tool for analyzing multi-mass tethered systems. The new model of the 3-mass 
tethered constellation has been used to perform simulations which show the capa- 
bilities and the limitations of the tether system with respect to the reproduc- 
tion of predetermined acceleration profiles on board the middle platform. Si- 
nusoidal vertical acceleration variations of the desired frequency and amplitude 
can be obtained at the expense of a smaller horizontal acceleration component. 
This component deflects the total acceleration from the local vertical during 
the sinusoidal motion of the middle platform. 

The ”g-tuning" capability has been also demonstrated. An open loop control 
law that provides a satisfactory dynamic response has been devised. The middle 
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mass is moved from the system C.M. to the desired location along the tether 
according to a modified hyperbolic tangent control law. Time constants for the 
control law have been derived in order to have a satisfactory transient re- 
sponse . 


3.0 PROBLEMS ENCOUNTERED DURING THE REPORTING PERIOD 


None 


4.0 ACTIVITY PLANNED EOR THE NEXT REPORTING PERIOD 

During the next reporting period further studies will be performed on the 
"g-tuning." We will investigate the case of displacing the middle mass from a 
position near the Space Station to the zero-g point in order to achieve a zero 
acceleration condition. We will also investigate the effect, upon the overall 
dynamics of the constellation, of a fast displacement of the middle mass from 
one tether tip to the other . 



